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Abstract. We investigate the nature of the phase transition of the ferromagnetic Potts model 
with invisible states. The ferromagnetic Potts model with invisible states can be regarded 
as straightforward extension of the standard ferromagnetic Potts model. The invisible states 
contribute the entropy, however they do not affect the internal energy. They also do not change 
the symmetry which breaks at the transition temperature. The invisible states stimulate a 
first-order phase transition. We confirm that the first-order phase transition with spontaneous 
g-fold symmetry breaking for q = 2,3, and 4 takes place even on two-dimensional lattice by 
Monte Carlo simulation. We also find that the transition temperature decreases and the latent 
heat increases as the number of invisible states increases. 



1. Introduction 

In statistical physics, it has been a central issue to clarify a relationship between the symmetry 
which breaks at the transition point and the nature of the phase transition such as the order 
of the phase transition and the universality class. In order to consider such a relationship, the 
standard Potts model is often adopted[TJ [2]. In the standard Potts model, spins have q values 
and its interaction is represented by Kronecker's delta. Suppose we consider only two-body 
interaction, the Hamiltonian of the standard Potts model is given as 




(1) 



where denotes the pair of the interacted sites. Note that the two-state Potts model is 

equivalent to the Ising model. Then the (/-state Potts model can be regarded as the standard 
generalization of the Ising model. In the standard (/-state ferromagnetic (Jjj > 0) Potts model 
on the two dimensional lattice, second-order and first-order phase transition with spontaneous 
breaking of (/-fold symmetry occur when q < 4 and otherwise, respectively. The standard 
ferromagnetic Potts model has sometimes succeeded to explain the nature of phase transition in 
complicated theoretical models and real materials [2j. For example, the phase transition of the 
absorption of 3 He on graphite and the orbital order of transition metal oxides such as LaMn03 
can be explained by three-state ferromagnetic Potts model [3J [H [S]. 

Recently, there have been found novel examples that the nature of the phase transition is 
not consistent with that of the standard (/-state Potts model^ [5]. For instance, a first- 
order phase transition with spontaneous threefold symmetry breaking takes place in some two- 
dimensional frustrated systems with a number of competed short-range interactions. This phase 
transition seems to be a counterexample of the concept of the standard ferromagnetic (/-state 
Potts model, since the standard ferromagnetic three-state Potts model on two-dimensional lattice 
has a second-order phase transition as mentioned above. Since mechanism of such a phase 
transition is not clear, from a viewpoint of statistical physics, it is important to understand why 
a first-order phase transition with threefold symmetry breaking occurs even on two-dimensional 
lattice. 

Purpose of our study is to construct a simple model which can explain such a nature of phase 
transition by modifying the standard Potts model and to investigate a microscopic mechanism 
of changing the order of phase transition. We found that the order of phase transition of three- 
state Potts model changes by adding invisible states[9]. In this paper, we study the detail of the 
phase transition of Potts model with invisible states and consider the effect of invisible states. 
The organization of this paper is as follows. In Section 2, we introduce the Potts model with 
invisible states and review some properties of this model. In Section 3, we show the mean-field 
analysis and result obtained by Monte Carlo simulation. In Section 4, we summarize our study 
and show future perspective. 

2. Model 

We consider the following Hamiltonian which is called the ((/,r)-state Potts model [91 110] : 

Q 

H = - J ^2 s si,sj ^2^s,, a S SjtCt , Si = 1,- ■ ■ ,q,q + 1,-- ■ ,q + r, (2) 

(i,j) Q=l 

where denotes the pair of the nearest neighbor sites on square lattice. Note that spins 

can have discrete values from 1 to q + r in this model, while the spins in the standard (/-state 
ferromagnetic Potts model can have a value from 1 to q. We call the states where 1 < Sj < q 
"colored states" and the states where q + l<8i<q + r "invisible states" [9]. In this paper, we 
only consider the case for ferromagnetic coupling J > 0. If and only if 1 < Sj = Sj < q, the 
interaction J works and it is an energetically favorable state. 

First we consider two spin system. Figure [1] (a) shows the ground states for the (3,2)-state 
Potts model. Figure Q] (b) and (c) denote excited states for the (3,2)-state Potts model. States 
depicted in Fig. [1] (a) and (b) appear in the standard ferromagnetic (/-state Potts model, whereas 
states depicted in Fig. [1] (c) appear only when the invisible states are introduced. Figure Q] (d) 
shows the energy level structure of two spin system for the ((/,r)-state Potts model. The number 
of ground states is q and that of the excited states is q 2 — q + 2qr + r 2 . Contribution of the 
invisible states in excited states is 2qr + r 2 . This energy level structure indicates that (/-fold 
symmetry breaks at the transition point. 
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Figure 1. (a)-(c): The (3,2)-state Potts model. The squares and dotted circles indicate colored 
state and invisible states, respectively, (a) Ground states of two spin system, (b) Excited states 
of two spin system. Both (a) and (b) also appear in the case of the standard ferromagnetic 
(/-state Potts model, (c) Excited states of two spin system due to invisible states, (d) Energy 
level structure of two spin system for the (g,r)-state Potts model. The number of ground states 
and excited states are q and q 2 — q + 2qr + r 2 , respectively. The dotted boxes in (c) and (d) 
denote the contribution of invisible states. 



By comparing the partition function, we can obtain equivalent Hamiltonian: 

where we rename the label of invisible state from q + 1 < Si < q + r to di = 0. The second 
term of Eq. ([3|) represents the chemical potential of invisible state. The temperature-dependent 
chemical potential comes from entropy effect of invisible states. Thus, the invisible state does 
not contribute the internal energy but affect the entropy. 

3. Result 

In the first half of this section, we consider the Bragg- Williams approximation. In the last half 
of this section, we study a nature of phase transition by Monte Carlo simulation. 

3.1. Mean-field analysis 

First we consider the order of the phase transition of the Potts model with invisible states by 
Bragg- Williams approximation. Let x a be the fraction of a-th state (0 < a < q) with condition 
Y^a=o x a = 1. To analyzing the phase transition with spontaneous breaking of g-fold symmetry, 
we assume x a as follows: 

x = t, xi = -(l-t)[l + (g-l) S ], x Q = -(l-t)(l-a) {2<a<q), (4) 

where < s, t < 1. Let z be the number of the nearest neighbor sites. The internal energy and 
the entropy are given by 

E BW ( s ,t) = - Z -±Y j xl- x T log r, 5 BW (a, t) = - J>„ log x a , (5) 
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Figure 2. The transition temperature T BW and the latent heat Ai? BW as functions of the 
number of invisible states r for z = 4 obtained by the Bragg- Williams approximation. 



respectively. Then the free energy is given by 
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(6) 



(7) 



From this free energy, we can obtain the latent heat and transition temperature in the cases of 
q = 2, 3, and 4. Figure [2] shows r-dependency of the latent heat and transition temperature for 
z = 4. From Fig. [21 we find a second-order phase transition takes place only for (2,1), (2,2), 
and (2,3)-state Potts model. Note that the transition temperature and the latent heat for r = 
(i.e. standard Potts model) are given as 



T C BW ( g = 2,r = 0) 
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(for g > 3), (8) 



AE BW (q,r = 0) 
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respectively [IT] . As the number of invisible states r increases, the latent heat AE BW increases 
and the transition point T BW decreases. From this fact, it can be considered that the invisible 
state drives a first-order phase transition. 



3.2. Monte Carlo simulation 

To investigate the effect of invisible states on two dimensional lattice, we study equilibrium 
properties of the (g,r)-state Potts model on square lattice by Monte Carlo simulation. Figure 
[3] shows typical snapshots of the (3,27)-state Potts model for several temperatures. The color 
points and white points indicate the colored states and the invisible states, respectively. The 
red, green, and blue points depict cij = 1,2, and 3, respectively. At zero temperature, all spins 
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Figure 3. Typical snapshots for the (3,27)-state ferromagnetic Potts model. As the temperature 
decreases, the density of invisible states decreases and the ordering domain grows. 



take the same colored state. As the temperature increases, size of the ordered domains become 
small. Above the transition temperature, small colored islands are floating in an invisible sea. 

Figure H] (a) shows the specific heat as a function of temperature. In order to obtain the data 
with high accuracy, we adopt reweighting method. The maximum value of specific heat increases 
almost linear dependence with system size. To clarify whether the first-order phase transition 
occurs or not, next we calculate the energy distribution and implement finite size scaling. Figure 
H] (b) shows the energy histogram at the temperature where the specific heat has the maximum 
value. Energy histogram in Fig. H] (b) is a bimodal distribution which is a characteristic behavior 
of the first-order phase transition. As the system size increases, the peaks become high. 

Next we implement the finite size scaling which is established by Challa et al. \12\. The finite 
size scaling for a first-order phase transition in d-dimensional system is treated as: 




T L=48(T=0.585280) L=96(T=0.585174) 
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Figure 4. (a) Temperature dependency of the specific heat for the (3,27)-state Potts model, 
(b) Energy histogram for the (3,27)-state Potts model. 



Table 1. Transition temperature and latent heat for several parameter set of (g,r)-state Potts 
model. For r = case, T c and AE were obtained exactly [2]. 
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Symmetry 


(2,0) 


1.13459 





2-fold 


(2,30) 


0.57837(1) 


1.02(2) 


2-fold 


(2,32) 


0.56857(1) 


1.23(2) 


2-fold 
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3-fold 
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3-fold 
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3-fold 
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0.61683(1) 


0.68(2) 


4-fold 
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4-fold 



where T c (L) denotes the temperature where the specific heat takes the maximum value C max (L) 
in finite size L d system. The values T c (oo) and AE denote the transition temperature and the 
latent heat for the thermodynamic limit, respectively. We fit T c (L) as a function of L~ 2 as 
shown in Fig.[5ja). From the intercept of fitting curve, we obtain T c (oo) = 0.58513(1). Next we 
fit Cmax {L) as a function of L 2 . We estimate the latent heat AE = 1.05(2) from the gradient 
of the fitting curve and the obtained transition temperature T c (oo) as shown in Fig. [5£b). 



0.5856 



0.5854 



0.5852 



0.585 




8000 



6000 



4000 



2000 



0.0005 

, -2 



0.001 




10000 



(a) 



(b) 



Figure 5. (a) T c (L) as a function of L 2 . Intercept of the fitting curve indicates the estimated 
transition temperature for the thermodynamic limit, (b) C max (L) as a function of L 2 . 

By the same analysis, we calculate the transition temperature and the latent heat for several 
parameter set of (q,r) (Table [1]). As the number of invisible states increases, the transition 
temperature decreases and latent heat increases. From Table [IJ we conclude the invisible states 
stimulate a first-order phase transition, which is consistent with the result obtained by Bragg- 
Williams approximation. 



4. Conclusion and Future Perspective 

In this paper, we study the nature of the phase transition of the Potts model with invisible 
states. The Potts model with invisible states is regarded as the straightforward extension of 
the standard ferromagnetic Potts model, since this model is constructed by just adding the 



invisible states. The invisible state does not affect the internal energy but contributes to the 
entropy. It does not change the symmetry which breaks at the transition temperature. From 
the result obtained by Bragg- Williams approximation and Monte Carlo simulation, the invisible 
state drives a first-order phase transition, since the latent heat increases as the number of the 
invisible states increases. We believe that this model is a fundamental model and it enables us 
to understand why the first-order phase transition with spontaneous 2,3, and 4-fold symmetry 
breaking can occur even on two-dimensional lattice. 

Recently, the Potts model has been adopted for analysis of problems in information science 
and technology such as clustering problem[T31 E]. It has been suggested that the effects of 
thermal and quantum fluctuation are efficient for optimization problems. We expect that the 
fluctuation which comes from the invisible states is efficient for such problems. 
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